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Abstract 

bX)| We examine the Unruh-DeWitt particle detector coupled to a scalar field in an 

arbitrary Hadamard state in four-dimensional curved spacetime. Using smooth 
' switching functions to turn on and off the interaction, we obtain a regulator- 

^ ^ , free integral formula for the total excitation probability, and we show that an 

instantaneous transition rate can be recovered in a suitable limit. Previous results 

■ in Minkowski space are recovered as a special case. As applications, we consider 
. an inertial detector in the Rindler vacuum and a detector at rest in a static 

Newtonian gravitational field. Gravitational corrections to decay rates in atomic 

■ physics laboratory experiments on the surface of the Earth are estimated to be 
O ' suppressed by 42 orders of magnitude. 

> '■ 

^ : 1 Introduction 

The Unruh-DeWitt model for a particle detector PQ [2] is an important tool for probing 
the physics of quantum fields wherever noninertial observers or curved backgrounds are 
present. In such cases there is often no distinguished notion of a "particle," analogous to 
the plane-wave modes in Minkowski space, but an operational meaning can be attached 
to the concept by analysing the transitions induced among the energy levels of a detector 
coupled to the field. Upwards or downwards transitions can then be interpreted as due 
to absorption or emission of field quanta, or particles. The best-known applications 
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of this procedure are those for which the spectrum of transitions is thermal, which is 
the case for uniformly accelerated detectors in Minkowski space [T], inertial detectors 
in de Sitter space [3j, and detectors at rest in the exterior Schwarzschild black hole 
spacetime 

In first-order perturbation theory, the transition probability of the Unruh-DeWitt 
detector is proportional to a quantity known as the response function, which involves 
integrating the Wightman distribution of the quantum field over the worldline of the 
detector. When the quantum state of the field is sufficiently regular and the detector is 
switched on and off smoothly, the response function is well defined [5], and the physical 
interpretation is that the response function is then proportional to the probability of a 
transition to have occurred by a time at which all interaction has already ceased. If, 
however, one wishes to address the probability of a transition to have occurred by a 
time at which the interaction is still ongoing, the response function is no longer well 
defined because the switching function then has a sharp cut-off at a singularity of the 
Wightman distribution. In special cases in which the trajectory is stationary, the vacuum 
state is invariant under the Killing vector generating this stationary motion and the 
detector has been switched on in the infinite past [H [21 El S] , the issue can be bypassed 
by formally integrating over the whole trajectory and factoring out the infinite total 
proper time, because by stationarity the transition rate can then be argued to be time- 
independent. But in a general setting this is not possible, and seemingly inconspicuous 
regularisations of the Wightman distribution can lead to unphysical results, even for 
uniformly accelerated motion in Minkowski space [6l[7j. 

A way to address this problem is to regard the sharp detector switch-off as a limit 
of a family of smooth switch-offs and investigate how the results depend on the way in 
which the limit is taken. In [8] this issue was investigated for a massless scalar field in 
four-dimensional Minkowski space, with the quantum field in the Minkowski vacuum. 
The response function with a smooth switching function was written in a form in which 
the integrand is no longer a distribution but a genuine function, and it was shown 
that a well-defined notion of a transition rate emerges when the switching time scale 
is small compared with the total duration of the coupling. It was also shown that in 
the appropriate limits this transition rate coincides with that obtained by regularising 
the sharply switched-off detector by a nonzero spatial size P, [7] . The key point is that 
when the Wightman distribution under the integral is represented by an ze-regularised 
function, the regulator limit e — > and the limit to sharp switching do not in general 
commute and the first must be taken before the second. 

The aim of this paper is to extend these results to a more general setting. For this, 
we will start in section [2] with a review of the Unruh-DeWitt detector, with special 
attention to the procedure introduced in [S] that allows limits of switching functions 
to be considered. In section [S] the results of [H] are generalised to a situation in which 
Minkowski space is replaced by an arbitrary four- dimensional globally hyperbolic space- 
time, the Minkowski vacuum state by an arbitrary Hadamard state and the massless 
scalar field by a scalar field with arbitrary mass and curvature coupling. We shall in 
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particular obtain a simple and manifestly well-defined expression for the difference in 
the response of detectors that have the same switching function but move in different 
quantum states of the field, on different trajectories or even in different spacetimes. The 
limit of sharp switching is discussed in section HI In sections [5] and [6] we use these results 
to obtain the detector transition rate in two examples of interest: an inertial detector 
in the Rindler vacuum in Minkowski space, and a detector at rest in a static Newtonian 
gravitational field. The results are summarised and discussed in section [71 Certain 
technical properties of the detector response in the Rindler vacuum are established in 
the Appendix. 

Throughout this paper we will assume a Lorentzian metric of signature ( — 
using the (+++) sign convention of Misner, Thorne and Wheeler [9]. We use units in 
which c = h = 1, while keeping G = Ip ^ 1. Spacetime points are denoted by sans-serif 
letters. The symbol 0{x) denotes a quantity for which 0{x)/x is bounded as a; — > 0. 
0(1) denotes a quantity that is bounded in the limit under consideration. 

2 Particle detectors and their regularisation 

We consider a detector consisting of an idealised atom with two energy levels, \0)d 
and \ l)d, with associated energy eigenvalues and u. The detector is following a timelike 
C°° trajectory x(r), parametrised by its proper time r, in a four- dimensional Lorentzian 
globally hyperbolic C°° manifold M. The coupling of the detector to a real scalar field 
of mass m and curvature coupling ^ is given by the interaction Hamiltonian 



where c is a coupling constant, /i(r) is the detector's monopole moment operator and 
x(t) is a smooth non-negative function of compact support, x is called the switching 
function: the interaction takes place only when x is nonvanishing, and because x 
compact support the interaction has a finite duration. If x(r) is not defined for all r G M, 
we assume the support of x{^) to be in the open interval in which x(r) is defined. 

We take the initial state of the joint system before the interaction to be |\I') ® \0)d, 
where the field state is an arbitrary Hadamard state [TUl |TT]. We are interested in 
the probability for the detector to be observed at state \l)d after the interaction has 
been switched off. Treating the coupling constant c as a small parameter, working to 
first order in perturbation theory in c, and summing over the unobserved final state of 
the field, this probability reads 0, |T21 [13] 




(2.1) 




(2.2) 



where the response function F{uj) is given by 



F{iu) = dr' dr' 



e-'-^^'-"h{r')x{r")W{T',r") 



(2.3) 



J — oo J — oo 
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and the distributional correlation function W{t', t") is the pull-back of the Wightman 
distribution W(x,x') := (\E'|0(x)0(x')|^) to the detector worldlineQ 

W{t',t") ■.= W{x{t'),x{t")) . (2.4) 

The response function thus encodes the properties that depend on the state |\E') and the 
detector trajectory, while the prefactor in fl2.2l) is a constant that only depends on the 
detector's internal properties. We shall from now on suppress the prefactor and refer to 
the response function simply as the probability. 

To summarise, (12.31) gives an unambiguous answer to the question "What is the 
probability of the detector being observed in the state \l)d after the interaction has 
ceased?" 

The meaning of the distributional correlation function under the integral in (12. 3p is 
somewhat subtle. Recall that the Wightman distribution W{x, x') in a Hadamard state 
can be represented by a family of functions fiUl ITT] 



WJx,x') 



(2vr) 



^ + t;(x,x')ln(a,(x,xO)+//(x,x') 



(TJX, X' 



(2.5) 



where e is a positive parameter, a{x, x') is the squared geodesic distance between x 
and x', a^(x,x') := a(x,x') + 2ie[T{x) — T(x')] + and T is any globally-defined future- 
increasing function. The logarithm denotes the branch that is real-valued on the 
positive real axis and has the cut on the negative real axis. A(x, x') is the Van Vleck 
determinant, which is smooth for sufficiently near-by x and x', the function f (x, x') is a 
polynomial in a{x,x'), and the function H{x,x') can be chosen for arbitrarily large 
n by taking the degree of the polynomial f(x, x') sufficiently high. The ie-prescription 
in (12.51) defines the singular part of W{x,x'): the action of the Wightman distribution is 
obtained by integrating We(x,x') against test functions and taking the limit e — 0, and 
this limit can be shown to be independent of the choice of the global time function T. 
Now, the distributional correlation function W{t,t') (12.41) is the pull-back of W{x,x') 
to the detector's worldline, which is a C°° submanifold. It follows that the action 
of the distribution W{t,t') is obtained by pulling back the function W^{x,x') to the 
function W^{t, t'), integrating W^{t, t') against test functions and taking the limit e — 
[5l [HI [15] . Formula (12. Sp must thus be understood as 

/oo /»oo 
dr' / dr" e~^-(-'--") x{'r')x{r") W,{r\ r") , (2.6) 
^oo J^oo 

where the integrand is now an ordinary function and the singular part of W{t, t') has 
been encoded in the ie prescription. As W^{x,x') = We{x',x), where the oveline denotes 



""^We denote both the spacetime Wightman distribution Vl^(x, x') and its pull-back W{t,t') to the 
detector worldline by the same letter, writing out the arguments explicitly in places where ambiguity 
could arise. 
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complex conjugation, we have W^{t',t") = W^{t",t'), and it follows that (12. 6p can be 
written in the equivalent form [6, _8J 

/oo noo 
duxiu) / dsxiu-s)e-''^'W,{u,u-s). (2.7) 
oo Jo 

Although formulas (12. 6p and (12. 7p are suitable for computing the detector's response, 
these formulas do not display a clear separation between those properties of the response 
that depend on on the trajectory and the quantum state and those properties that only 
depend on the choice of the switching function. Neither do these formulas exhibit how 
the response depends on the proper time along a given trajectory. Several authors 
[SI [ISl [13, [IE] have therefore addressed the question: "If the detector is turned on at 
proper time Tq and read at proper time r, while the interaction is still on, what is the 
probability that the transition has taken place?" If issues of regularisation could be 
ignored, this would amount to adopting in (12. 3p the switching function 

x(r') = e(r'-ro)e(r-r'), (2.8) 

where G is the Heaviside function. The transition probability then becomes a function 
of the reading time r and can be written as 

Fr{uj) = 2Rei du dse-''^'W{u,u- s) , (2.9) 

Jto Jo 

and we can define the instantaneous transition rate as its derivative with respect to r, 

/•At 

Friuj)=2Re dse-''^'W{T,T- s), (2.10) 
Jo 

where At := r — tq. Ft-{uj) thus represents the number of transitions per unit time 
in an ensemble of identical detectors. It is this instantaneous transition rate that one 
expects to have the r-independent Planckian spectrum in the Unruh effect in Minkowski 
spacetime and in its generalisations to curved spacetimes, once tq is taken to — oo to 
avoid transient effects. A comprehensive recent review of the Unruh effect can be found 
in pL9j. 

We note in passing that the transition rate (I2.10p is not directly related to transition 
rates that could be measured with a single ensemble of detectors. Given an ensemble 
of identical detectors on a given trajectory, -FV(^) gives the fraction of detectors that 
have undergone a transition when observed at time r, but as an observation alters the 
dynamics of the system, Ft-{uj) no longer has this interpretation after a first observation 
has been made. To measure ^^-(a;), one therefore needs a set of identical ensembles, 
such that each ensemble is used to measure F^-^u) at just a single value of r. Note in 
particular that Ft{^) may well be negative at some values of r [201 EI]- -^r(t^) may thus 
be difficult to measure operationally, but it is nevertheless of interest as a nonstationary 
generalisation of the transition rate that naturally arises in stationary situations. 
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Returning to formulas fl2.9p and (12.101) . the difficulty with them as written is that 
the 'switching function' fl2.8p is not smooth. We are no longer guaranteed that replacing 
W{t, t') by W^{t, t') in (12. 9p and fl2.10p and taking the hmit e — > would give a result 
that is independent of the choice of the global time function in (12. 5p . Case studies have 
shown that the result depends on the choice of the time function for Minkowski vacuum 
in Minkowski space [6l [71 [20l [21] and the Euclidean vacuum in de Sitter space [201 [21] , 
and the methods of Appendix A of [7] can be adapted to show that the same holds for 
arbitrary Hadamard states in an arbitrary spacetime. Formula (I2.10p does therefore not 
provide a well-defined notion of an instantaneous transition rate. 

One way to address this problem was introduced in [6] and further developed in [711201 
2T] . The idea is to replace the correlation function W{t,t') in (I2.10p by a correlation 
function in which the field operator has been smeared over a spacelike hypersurface 
orthogonal to the trajectory. The weight function in the smearing is characterised by 
a positive length parameter e, which acts gulator and corresponds physically to 

the spatial size of the detector in its instantaneous rest frame. At the end the pointlike 
detector limit e ^ is taken. This scheme does not rely on the choice of a time function 
to regularise the Wightman distribution, and in Minkowski space the introduction of 
the spatial hypersurfaces is straightforward [B] and there are partial results regarding 
independence of the choice of the weight function [7]. An implementation of the scheme 
in de Sitter space was given in [201 [211 [22] . However, the spacelike surfaces introduced in 
Minkowski space in [6j are not easily generalisable to spacetimes without a high degree 
of symmetry, and it would seem desirable to attach a meaning to the instantaneous 
transition rate within the framework of the conventional regularisation of the Wightman 
distribution (12. 5p . 

A way that stays fully within the conventional regularisation of (12.51) . (12.61) and (12.70 
was introduced in |8] in the special case of Minkowski spacetime, massless scalar field 
and the Minkowski vacuum state. Adopting a Lorentz frame with global Minkowski 
coordinates (t, x) and choosing t as the global time function, formula (12. 7p for the 
transition probability becomes 



OO fOO 
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where (Ax)^ is the squared geodesic distance between x{u) and x{u — s) and At 
t{u) — t{u — s). The limit e — > can be computed explicitly, with the result [H] 



+ — / dux{u) dsx{u-s)[-^ + -] . (2.12) 
27r^ Jo V (Ax) 



^(^) = / du[x(u)]2+^/ ^/ <\ux{u)[x{u) -xiu- s)] 



When the switching function x equals 1 over an interval of length Ar, and the switch-on 
and switch-off each take place within an interval of length 5 with a profile that scales 
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with 5 but whose shape is otherwise fixed, the leading behaviour of the transition rate 
(defined as the derivative of F{u) (12. 121) with respect to At) at 5 is 

Friu^) = + r ( + \) +^ + 0{S) , (2.13) 

where now (Ax)^ is the squared geodesic distance between x(r) and x(r — s). In the hmit 
(5 — >■ 0, the transition rate (I2.13P agrees with that obtained from spatial smearing in [7], 
and it reproduces the expected Planckian spectrum when the trajectory is uniformly 
linearly accelerated and Ar oo. Further properties of this transition rate are discussed 
in [ZIIH]. 

In this paper we generalise the Minkowski vacuum results fl2.12p and (I2.13P to a 
general Hadamard vacuum state in four- dimensional spacetime, for a field with arbitrary 
values of the mass and the curvature coupling. We shall show that most of the arguments 
in [8] carry over to this situation, and we shall find the expressions that generalise (12.121) 
and (12.131) . These expressions will then be applied to two examples. 

3 Regulator-free response function in a general 
Hadamard state 

In this section we obtain a regulator-free expression for the response function F{u}) by 
computing explicitly the limit e — >^ in (12. 7p . Following the procedure used in [S], we 
split the s-integral into the subintervals (0,//) and {rj, oo), with r] = y/e, estimate the 
integrand in each subinterval and finally combine the results. 
We shall make use of the small s expansions 



a 


= -s'-Lah' + 0{s'), 


(3.1a) 


A 


= l + 0{s'), 


(3.1b) 


V 


= m'+(^^-^^R + 0{s'), 


(3.1c) 


AT 


= Ts-'-^ + 0{s'), 


(3.1d) 



where the Ricci scalar R, the squared (covariant) acceleration and T are evaluated 
at the point x{u) and the dots indicate proper time derivatives. 

3.1 Subinterval s G (r/, oo) 

Consider in (12. 7p the subinterval s G {rj, oo), and let Wq denote the pointwise limit of 
as e — 0. Replacing by Wq creates under the w-integral an error that equals x{u) 
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times the quantity 

I'OO 

2Re / dsx{u - s)e~''^'[W,{u,u - s) -Wq{u,u - s)] 

Jri 

= ^ Re ds x{u - s) e''-^ | A (^1 _ + ^ [in(a,) - ln(a)] | , (3.2) 

where the functions A, v, a and are each evaluated at the pair (x, x') = (x(u), x(m — s)) 
and ln((T) := hm£^o+ ln(o"g). We shall show that this error term does not contribute to 
(12.71) after the limit e — > is taken. 

Consider first in (13.21) the contribution from o"~^ and a~^. We split this contribution 
into its odd and even parts in uj. The part that is odd in uo can be written as 

1 2eAT 

Jn . f f^V 4e2(AT)^ 



27r' 



1 + ^1 + 



a J 0"^ 



where AT := T(x(m)) — T(x('U — s)). As the switching function makes the upper limit of 
the s-integral finite, it follows from (13. lap and (]3.1dl) that the quantities (AT)^/o" and 
e/o" are bounded by constants that are independent of rj. We can therefore write 

as 

1 /"°° , / N A 1 /9 e sin(co's)AT r ^, .1 , 
— dsxiu-s)A'/^ V [l+0(e)], (3.4) 

^ J r) ^ 

where the 0(e) estimate holds uniformly in s. As the functions x ^"^d A are 0(1) at 
small s, the integrand in (13. 4p is bounded by a constant times e and (13. 4p is thus of 
order 0{eri~^) = 0{ri). Similarly, the part that is even in uj can be written as 

6^ 4(AT)^ 
^2 1 + 



, , ds xiu - s) cos(cus) Ai/2 a a (3 5) 



(J I 0"^ 



and similar estimates show that the integrand in (13.51) is bounded by a constant 
times e^s~^. The expression (13. 5p is hence of order 0(e^?7~^) = 0{r]). 

Consider then in (13. 2p the contribution from the logarithmic terms. Keeping track 
of the branches of the logarithms, we can write this contribution as 

Re / dsx[u - s) e-'^'v In 1 + + - • (3-6) 



27r^ jj? V (J ^ 

It follows from (13. lap and ( 13.1dp that eAT/a is bounded by r] times a constant and e^/cr 
is bounded by e times a constant. The logarithm is hence of order 0{t]) uniformly in s, 
and (13. 6p is of order 0(77). 

As X has compact support, all the estimates above hold uniformly in u under the 
M-integral in (12. 7p . In the subinterval s G [i], 00) in (12.70 . can therefore be replaced 
by Wq without error. 
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3.2 Subinterval s G (0, ij) 

We now turn to the subinterval s G [0,1]) in fl2.7p . The singularity of W at s = 
implies that we cannot directly replace by Wo, and we will need to examine the 
small s behaviour of more closely. 

We observe first that the term H{x, x') in (12.51) clearly gives a vanishing contribution 

to (122D. 

Consider then the logarithmic term in (12. 5p . Suppressing for the moment the fac- 
tor the integral over u and the limit e ^ 0, the contribution to (12. 7p reads 

^ \0{s')]\n{a,). (3.7) 



27r2 



Re / ds xiu — s)e~ 



R 



The imaginary part of the logarithm is bounded and its contribution in (13.70 is therefore 
of order 0{ri). To estimate the real part of the logarithm, we write s = ex, with 
< X < l/r], and use the expansions (13. ip to obtain 



'l-x'f + 4x'f' 



[l + 0{v)] 



where the 0{ri) term holds uniformly in x. Hence 



2Reln((Te) = 41ne + ln 



X 



2\2 



+ Ax'^f'^ 



+ 0{v) 



(3.8) 



(3.9) 



where again the 0{ri) term holds uniformly in x. The contribution in (13.70 is therefore of 
order 0(77 In r/). As this estimate holds uniformly in u, by virtue of the compact support 
of the logarithmic term does thus not contribute in (12.70 . 

Finally, consider the term in (12. 5p . From (13.80 we see that s"^ is bounded, 
and hence ds s"^ = 0(7]). It follows from ( ]3.1bO that we may replace A^/^ by 1, 
and we may similarly replace the factor x(m — s)e~*'^'' by (1 — iujs)x{u) — sx{u). What 
remains is to analyse the small 77 behaviour of the expression 



27r2 



Re 



''^^ (1 - i^s)x - sx 



(3.10) 



where x ^"^^ X ^"^^ evaluated at u. In the special case in which the spacetime is Minkowski 
space and the global time function T is the Minkowski time coordinate in a given Lorentz 
frame, this analysis was carried out in [8], and the techniques used therein generalise 
to (13. 100 in a straightforward way. Splitting /< into its even and odd parts in uj as 



J< = I'^^"^ -\- J°^'^, and writing s = ex with < x < l/?7, we fincj^ 



27r2 



1/7, 



1 — x^) dx 



'1 - X2)' +4x2T2 



X 



rodd 



^XT 



x^ dx 



[l-x^Y + Ax^T^ 



XX + 



AxTTx^ 



X 



2)" + 4x2T2 



+ 0{ii), 

(3.11a) 
(3.11b) 



^Our formula (j3.11bp corrects a typographical error in formula (3.4b) of [5]. 
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In (13. Hap the integral of the first term is elementary, and multiplying the second and 
third term by x yields a total w-derivative that can be taken outside the integral. The 
result is 



.even ^ ^ _ ^± <l ' 

^ < 27r2r7 47r2 du J, (1 - x^)^ + Ax^f^ 
The integral in (13.11bp is elementary, and multiplying the result by x we obtain 

xit' = + 0{v) . 



(3.12) 



(3.13) 



All these estimates hold uniformly in u, owing to the compact support of x- The only 
terms that contribute in the subinterval s G (0, t]) in (12.71) are therefore the explicitly- 
displayed terms in (13.121) and (I3.13p . 



3.3 Joining the subintervals 

Substituting the results of subsections 13. II and 13.21 in (12.71) . we find 



F(^) 



47r 



+ 21im 



du [xiu)f 
duxiu) 



dsxiu- s)e-''^'Wo{u,u- s) 



. (3.14) 



As X has compact support, the total derivative term in (I3.12p integrates to zero and has 
dropped out. What remains is to take the limit in (I3.14p . 

Following [8j, we take the term proportional to l/?7 under the s-integral, add and 
subtract under the s-integral the term x{u — s)/ (47r^s^) and group the terms in the form 



^j_Ju[x{u)f 
du xiu) 



+ 2 lim 



ds Re 



xiu-s)e-''^'Wo{u,u-s) + 



+ lirn 2 / duxiu) / dsxiu — s)Re 



+ 



1 



27r2./, s 



V 

^ ds 



e-'^'Woiu,u- s) + 
dux{u)[x{u) - x{u - s)] \ , 



4^2_g2 



(3.15) 



oo 



where in the last term the interchange of the m- integral and the s-integral is justified 
by absolute convergence of the double integral. The limit r/ — > can now be taken 
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by simply setting 77 = 0. In the last term the reason is that the ti-integral, when 
regarded as a function of s, has a Taylor expansion that starts with O(s^). In the 
term involving Wq, the reason is that the real part of e"*"^^ Wq{u, u — s) has the small s 
behaviour of — l/(47r^s^) plus an integrable function of s, by virtue of (12.51) and f l3.1al) . 
The final result for the response function is thus 

io 1 f°° ds f°° 

^(t^) = J du[xiu)f + ^ j dux{u)[x{u) -x{u- s)] 



/oo /»oo 
dux{u) I dsx{u — s)B.e 
00 Jo 



e-^'^'Waiu.u- s)+ ^ 



(3.16) 



In the special case of the Minkowski vacuum in Minkowski space, (13.161) duly reduces 
to the expression (12.121) found in [8j. 

The first two terms in (I3.16P depend only on the switching function x but neither 
on the quantum state, the spacetime or the trajectory. If we compare two detectors 
in different quantum states of the field, on different trajectories or even in different 
spacetimes, but having the same switching function, the difference of the responses is 
given by 



\u,u — s 



/oo /"OO 
dw x{u) I ds x{u - s) e~''^' \W^{u, u - s) - 
-00 Jq 

(3.17) 

where Wq and Wq are the pull-backs of the unregularised Wightman distributions in the 
two situations. The representation (12.51) of the Wightman distribution in a Hadamard 
state guarantees that the divergences in (I3.17P cancel and the integral is well defined. 
This is particularly convenient for numerical calculations. 



4 Sharp switching hmit 

In this section we discuss the response function (13.161) in the limit of sharp switch-on 
and switch-off. As in the case of Minkowski vacuum p], we shall isolate the divergence 
due to the sharp switching from a finite remainder and show that a well-defined notion 
of instantaneous transition rate can be defined in an appropriate limit. 

To control the switch-on and switch-off, we assume the switching function to have 
the fornix 

= ft, (41) 

where r, tq and 5 are parameters satisfying tq < r and < 5, and /ij, i = 1,2, are 
non-negative C°° functions satisfying hi{x) = for x < and hi{x) = 1 for 1 < x. This 
means that the detector is turned on smoothly during the interval (tq — 5, tq), with a 
profile determined by the function hi, it then remains turned on at constant coupling 

■^Our formula ()4.ip corrects a typographical error in the argument of hi in equation (4.1) of [H]. 
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strength for the time Ar := r — Tq, and it is finally turned off smoothly during the 
interval (r, r + 6), with a profile determined by the function /i2. The functions hi are 
regarded as fixed. We initially regard Tq and r as fixed but will eventually allow r to 
vary. 

The first term in (13.161) is equal to — (ct;/47r)(Ar + 6Ci), where Ci is a positive con- 
stant. The second term in (13.161) was analysed in [8], with the result that it only depends 
on Ar and S through the combination 6/ At and has at small 6/ At the asymptotic form 



1 [At 



Co 



(4.2) 



where C2 is a constant and the full expansion of the 0-term proceeds in positive powers 
of 5/ At. The last term in (I3.16P can be analysed by breaking the integrations into the 
various subintervals as in [S], with the result 



2 / du 

'TO Jo 



U — TO 



ds Re 



Wo{u,u- s) + 



+ 0(6). 



(4.3) 



The qualitatively new feature compared with P] is the logarithmic singularity in Wq, 
but the contribution from this singularity can be verified to be of order 0(5^ In 5) and 
hence subleading in (14. 3p . Collecting, we find 



UJ 



F(uj) = Ar 



U-TO 



+ 



27r2 



In 



+ 2 du 

J TO Jo 

At 



ds Re 



e-''''Wo{u,u- s) + 



1 



Att'^s^ 



+ C2 + 0{6) 



(4.4) 



By the smoothness of the spacetime and the trajectory, we may differentiate the 
trajectory-dependent contribution (14. 3 p with respect to r termwise, and the same holds 
for the trajectory-independent contributions by their explicit structure. We may there- 
fore take the r-derivative in (14. 4p termwise, with the result 



At 



ds Re 



e-^-Wo(r,r-s) + 



+ 



27r2Ar 



+ 0{6). (4.5) 



Equations (14. 4p and (14. 5 p are our main result. The transition probability (14. 4p di- 
verges as (5 — 0, but the divergence has been isolated into an explicit logarithmic term 
that is independent of the trajectory or the quantum state of the field. The r-derivative 
of the transition probability is given by (14. 5 p and remains finite as 5 — >■ 0. Equation (14.50 
provides a definition of what is meant by the detector's transition rate, without the need 
to introduce spatial profiles or other regulators. In the special case of a massless field in 
Minkowski spacetime, in the Minkowski vacuum, the limit 5 — > in (14. 5 p recovers the 
transition rate obtained in via a spatial profile regularisation. 
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We end the section with two comments. First, if we interpret the naive expression 
fl2.10p for the transition rate as 

/•Ar 

F,(cu) = hm2Re / ds e-''^' W,{t, r - s) (4.6) 
Jo 

and apply the methods of section [HI we find that (14.61) is equal to the 5 limit of 
our transition rate (14.51) plus an additional term proportional to T. The additional term 
vanishes if the pull-back of T to the trajectory is an affine function of r. The transition 
rate of a sharply-switched detector can thus be calculated equivalently from the 6 —>■ 
limit in (14. 5p . where the singularity in the Wightman distribution is cancelled by an 
explicit counterterm, or from the e — > limit in (14. 6p . provided the time time function 
used to regularise (14. 6 p is an affine function of r on the trajectory. 

Second, if we compare two detectors in different quantum states of the field, on 
different trajectories or even in different spacetimes, but having been in operation for 
the same length of proper time, we recover for the difference of the transition rates in 
the S ^ limit the formula 

AF,(cj) = 2 Re / ds e"^"^" [W^ir, r - s) - W^{t, t - s)] . (4.7) 
Jo 

The difference in the transition rates in two given situations can thus be written as a 
Fourier transform of a function that requires no regularisation. Formula (14. 7p is useful 
for both analytical and numerical calculations, especially in cases where the transition 
rate in one of the two situations is already known. In the next two sections we shall 
apply this formula to two such examples. 



5 Inertial detector in the Rindler vacuum 

In this section we consider a detector moving inertially through the Rindler wedge in 
Minkowski space, coupled to a massless scalar field in its Rindler vacuum state. This 
is the state in which the uniformly accelerated detectors associated with the Rindler 
wedge do not get excited. A naive apphcation of the equivalence principle could be 
argued to imply that as an accelerated detector moving through the 'unaccelerated' 
(Minkowski) vacuum state gets excited thermally, an unaccelerated detector moving 
through the 'accelerated' (Rindler) vacuum should also get excited thermally. Working 
in the limit 5 — 0, we shall show that this does not hold: the detector does have a 
nontrivial transition rate, but the rate is neither thermal nor constant in the detector's 
proper time, and it diverges as the detector approaches the Rindler horizon. 

Let {t, X, y, z) be a set of standard Minkowski coordinates in Minkowski space. We 
take the detector to move in the 'right- hand- side' Rindler wedge, x > denoted by R. 
The Rindler vacuum Wightman distribution iy^(x, x') in R reads [23] 

1"^ dv 

= W''{x,x') - / -^^W^'{xy{v)) , (5.1) 
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where W^{x,y) = [47r^(x — y)^] is the Minkowski vacuum Wightman distribution, 
the points x = {t,x,y,z) and x' = (t' , x' , y' , z') are in R and x"(t>) := (— t'coshf — 
x' sinhf , — x' coshf — t' sinhf , J/', z') . The difference of iy^(x, x') and W^''{x,x') consists 
thus of the integral of over an orbit of the associated Kilhng vector in the opposite 
Rindler wedge. Note that since the two Rindler wedges are spacehke separated, the 
difference term is a nonsingular function, and the distributional character of ^^'^(x, x') 
comes entirely from the first term on the right-hand side in (I5.ip . As we shall be using 
formula fl4.7p . we are suppressing the distributional issues in fl5.ll) . 
We take the trajectory of the detector to be 

x(r) = (r,X,0,0), (5.2) 

where X is a positive constant and r is the proper time. The trajectory stays in R during 
the proper time interval |r| < X. We must therefore consider the detector response in 
a finite proper time interval. 

We shall compute the transition rate from formula (14.71) . using the Minkowski vacuum 
as reference state. Unlike the infinite At case, in which the excitation rate is zero, a 
detector moving inertially through Minkowski vacuum over a finite proper time At does 
have transient excitations due to the switching. This transition rate has been found in 
several papers (see e.g. [IB]) and equals 

Fl'» = -£ + ^f^ + ^-S,(.Ar) (5.3) 

47r ZTT^Ar zvr^ 

where Si is the sine integral function. F^^^{uj) diverges as At (which is an artifact 
of omitting the 0{6) terms in (14.51) ) and approaches [to'/(27r)]0(— cu) for large At. 

Using (14. 7p with (15. ip . and substituting the trajectory (15. 2p with —X < tq < r < X, 
the transition rate in the Rindler vacuum at time r can be written as 

F^^'ico) = F^iico) + AF^co) , (5.4) 

where 



AFJuj) = —-—7 / ds cos(ti;s 



At 



dv 



1 



— [r + (r — s) coshf + Xsinhf]^ + [X + Xcoshf + (r — s) sinhf]^ 

(5.5) 

The u-integral in (15. 5p can be done by contour integration. Closing the contour in 
the upper half-plane, there are two infinite series of contributing poles, at respectively 
V = \n{X+T)-\n{X+T-s) + {2m+l)iTT andv = ln(X-r+s)-ln(X-r) + (2m+l)i7r with 
m = 0, 1, 2, . . ., and a single contributing pole at f = in. Summing over the residues, 
we find 



1 f^^ , cos(ti;s) 



AFr(u) = — 7 / ds 



27r2 7o s2 



1-:^ ^^^-^ + 



(5.6) 



14 



Note that the integrand in (15.61) remains finite as s — > 0, and the integral is well defined. 
Regarding tq fixed and r as variable, we see that AFt-{lj) tends to zero as r — *• Tq, but 
it diverges to — cxo as r — > X, which is the limit in which the trajectory approaches the 
Rindler horizon. We shall show in the Appendix that the asymptotic form of AFr{Lj) 
at r — > X is 



27r2X 









+ - In 




2 



In 



r 



0(1] 



(5.7) 



and the divergence is thus logarithmic in r. 

The response of inertial detectors in the Rindler vacuum was previously studied by 
Candelas and Sciama [21], but in a somewhat different framework. Candelas and Sciama 
investigate the whole family of trajectories (15.21) . and they compute the transition rate 
on each trajectory at the point r = — a^^, where a is a fixed positive constant. 
This means that the inertial trajectories are being compared along a Rindler trajectory 
of acceleration a. In the limit X —>■ oo, it is found that the response approaches the 
Minkowski vacuum value [ci;/(27r)]0(— tu). Candelas and Sciama's interpretation of this 
result is that "for the case of a charge moving inertially in Minkowski space-time through 
an accelerated vacuum the spectrum of field fiuctuations perceived by the charge is the 
same as if the vacuum were unaccelerated" (p4J, p. 1717). They view the limit of large 
X as a means to eliminate the transient effects due to the detector starting its inertial 
motion at r = 0. 

In our view these transient effects are already contained in the Minkowski vacuum 
part (15. 3p of the response, and the additional Rindler vacuum contribution (15. 6p shows 
that the inertial detector in the Rindler vacuum responds genuinely differently than in 
the Minkowski vacuum. In particular, the additional Rindler vacuum contribution (15.61) 
diverges as the trajectory approaches the Rindler horizon. This divergence was found 
previously by Davies and Ottewill [2Sj, both by a numerical evaluation of the response 
and by a comparison of the response with the expectation value of 0^. The divergence 
can indeed be expected on the grounds that in the Rindler vacuum the expectation 
values of both 0^ and the stress-energy tensor diverge on the horizon [231 ESI ES]. 



6 Detector at rest in Newtonian gravitational field 

In this section we shall find the response of an Unruh-DeWitt detector at rest in a static, 
Newtonian gravitational field: a static, asymptotically flat spacetime that satisfies the 
linearised Einstein equations with pressureless matter as source. The quantum field is 
assumed massless, but with arbitrary curvature coupling, and its vacuum state is taken 
to be the Boulware-like vacuum defined in terms of the global timelike Killing vector. 
We shall see that in this situation the detector's excitation rate is always zero, but the 
de-excitation rate (response for negative u) in general differs from that of an inertial 
detector in the Minkowski vacuum in Minkowski space, and the gravitational correction 
depends on the details of the mass distribution even when the detector is far from the 
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source. We also provide an order-of-magnitude estimate for this gravitational correction 
in atomic physics decay rates on the Earth's surface. 



6.1 General matter distribution 

The metric takes the form g^^, = 77^,^ + h^^, where t^^j, is the Minkowski metric and h^i, 
is the linearised correction. We use a system of Minkowski coordinates {t, x) in which 
rj^i, dx'^ dx" = — dt^ + dx^ and the correction components h^^, are independent of t. We 
assume /i^i^ to be small enough for validity of linearised Einstein's equations, and at 
large r := |x| we assume h^i, to have the asymptotically flat falloff 0(r^^). We further 
take hfj,i, to be in the Lorentz gauge, d^h^" = 0, where h^u := h^i, — ^rj^^h'^a- 

We need the Wightman distribution in this spacetime, in the Boulware-like vacuum 
that reduces to the Minkowski vacuum at the asymptotically fiat infinity. Working in 
perturbation theory to the order that is consistent with linearised Einstein's equations, 
this Wightman distribution must be the sum of the Minkowski vacuum contribution 
[47r^(Ax)^] and a correction W^^\x,x!) that is first-order in h^^, and dies off at infinity. 
Restricting the attention to a detector that remains at constant x and was switched on 
in the infinite past, equation (14. 7p shows that the transition rate reads 

F{uj) = e(-cu) + 2Re / dse-*"'P^W(t,x;t- s,x) . (6.1) 

2vr Jo 

The first-order Wightman distribution W^^"* in (16. ip depends on x and s but not on t, and 
from (I3.1cp we see that it has at s = an integrable logarithmic singularity proportional 
to the Ricci scalar. 

To find Vr(^) (x, x'), we first calculate the Feynman Green's function Gf(x, x') to first 
order in h^^, adapting the procedure that was introduced in [27] in the context of vacuum 
polarisation. We then find W'^^\x^x!) from the relation 

iG^(x, x') = W{x, x!) 0(t - f) + W{x!, x) e{t' - t) , (6.2) 

which is reliable order by order in perturbation theory as long as no new singularities 
turn up, and this will be seen to be the case for the first-order contributions. 
Consider the equation satisfied by G'i?(x, x'), 

[a - m^)] Gf{x, x') = <5(x, x') , (6.3) 

and expand g^^^, = -q^y -\- h^j, and Gf{x,x') = Gp\x,x') + Gp\x,x'), where Gp^ is the 
Minkowski vacuum Feynman propagator, 

GfM = . .r ..2~\ ^^-T ' (6-4) 



47r2 [-{t - t'f + |x 
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with its distributional part specified by the prescription e — 0+. Dropping second order 
terms and noting that G^p^ satisfies the zeroth-order equation, we obtain 

= [ Wd^ + e R^'^ (x)] op (x, x') , (6.5) 

where the last equality follows because in the Lorentz gauge hf^i, satisfies the Minkowski 
space wave equation 



-IGirGT, 



(6.6) 



and we are assuming static, pressureless matter. 

The key observation now is that because the spacetime is static and asymptotically 
fiat, Gp can be regarded as the analytic continuation of the unique Green's function on 
the positive definite section. As this holds order by order in perturbation theory, we can 
solve (16. 5p for G^ by using gP as the inverse of n^^\ with the result 



GPix,x') = J dxGPix,x) [h''ix)d^ + ^R^'\x)] GPix,x') 



167r4 



dx 



(t — t )^ + |x — x|^ + ie 



~ |2 



X 



[F\x)d? + ^R^'\k)] 



r^2 



{V - ty + |x'-x|^ + ie 



(6.7) 



The integral over t can be done by residues, using the time- independence of R^^^ and 
Defining gP{s, x) := GP{t, x; t — s, x), and writing X := |x — x| for short, we obtain 



gP{s,^) 



di 



s2 - 4 (X2 + ie] 



+ 



2/i°°(x) [3s^ + 4{X^ + te)] 
[s2-4(X2 + ze)]' 



. (6.8) 



As V^/iOO = -2i?(i), where := d^^ + dl^ + dl^, and as h^^ has the falloff O(r-i), 
we can integrate the term involving hP'^ by parts and take i?'-^^ as a common factor. 
Assuming s 7^ 0, and dropping terms that go to zero as e — > 0, we obtain 



gP{s,^) 



dii?«(i 



2e 



+ 



1 



2X [(s2 - Ait) - 4X2] x(s2 - Aie) 
1 



+ 



(s2 - Aie) 



3/2 



In 



Aie 



2X + Vs2 - Aie 



(6.9) 



Note that the 0{r ^) falloff of h^y implies that R^^^ has the falloff 0{r '^), and the 
integral in (16. 9p hence converges pointwise in s. As a consistency check, we also note 
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that if i?*^^^ vanishes in a neighbourhood of the point x, taking in (16.91) the hmit s ^ 
and e ^ yields a nonsingular expression that reproduces the vacuum polarisation 
(0^(x)) that was found for this class of spacetimes in \27\ . 

We now assume s > 0. Let Gp\{s,x.) and Gp\{s,x.) denote the contributions to 
G^p\s, x) from respectively the first term and the last two terms in the integral in (16.91) . 
Taking the limit e ^ in G^p\{^^^) is elementary, with the result 



X X"^ 
dXR{X)l- + — 



In 



\2X 



2X + S 



(6.10) 



where R{X) denotes the average of over a sphere of radius X about x and the 
dependence of R{X) on x is suppressed. In Gp\{s,:x.), splitting the integrand into 
partial fractions and taking the limit e — yields 



Gp\{s, x) 



2^ 



32n 



R{s/2) + i{HR){s/2) + i{HR){~s/2) 



where H stands for the Hilbert transform, 

{Hf){x) :=-F 

TT 



fjy) 

y-x 



(6.11) 



(6.12) 



with P denoting the principal value integral, and R is understood to vanish for 
negative argument. Recall now that the Hilbert transform can be written as 
Hf = —if+ + if-, where /+ and /_ are respectively the projections of / to 
the positive and negative frequency subspaces, f+{x) 



[x) := (27r 

-1/2 



-1/2 rO 



I. 



- (27r)-^/7;°e— /Mdo., 
e~^'^^ f {uj) duj , and / denotes the Fourier transform, /(u;) = 



(2n) ' e^^^f{x)dx [28]. As R vanishes for negative argument, it thus follows 
from flOT|) °that 



Gp\{s, x) 



2e 



167r 



R+is/2) + R+{-s/2) 



(6.13) 



Let now AFx(^^) '■= F(uj) + [(c(j/(27r)] 0(— cu) denote the correction to the Minkowski 
space transition rate. Using (16.11) and (16. 2p . we can write AFx(i^) in terms of Gp\ and 



G^p 2 as 



AFJu] 



2 Re / dse" 
'o 



iG^P,{s,^) + tG'^>,{s,^) 



(1) 



(6.14) 
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The contribution to AFx(co') from G''p\ equals 



2^- 


1 


Stt 




2e- 


1 


Svr 




2^- 


1 


Stt 






1 


47r 




2^-1 





Im / dse-*'^" /?+(s/2) + i?+(-s/2) 
Jo L 

/oo 
dse"^^" R+{s/2) 
-oo 

/oo 
dse-*^'^(s/2) 
-oo 

POO 

e(-cu)Im / dXe-2^'^^i?(X) 



ie(-.) |dxfl<')(x)2!H|££), 



(6.15) 



where we have used the definition of the positive frequency projection, changed the 
integration variable to X = s/2 and finally written R in terms of R^^\ To evaluate the 
contribution to AFx(co' 

by the absolute convergence of the double integral, and obtain 



from Gp\, we interchange the integrals over s and X, justified 



1 



oo 



— / dXR{X) 
^ Jo 



X Re 



oo 



dse" 



X^X^ 



In 



\2X 



2X + S 



(6.16) 



The integral over s in (16.161) may be interpreted as the integral of a complex analytic 
function along the positive real axis, with a contour deformation to the lower half-plane 
near the logarithmic singularity at s = 2X. For u > 0, the contour can be deformed to 
the negative imaginary axis, and the integral vanishes on taking the real part. For a; < 0, 
the contour can be deformed to that shown in Figure [H The contribution from the large 
arc vanishes when the arc is taken to infinity, and the contribution from the positive 
imaginary axis vanishes on taking the real part. The only nonvanishing contribution 
comes from the branch cut at s > 2X. Collecting, we find 



AF^,,{u;) = -&{-uj) rdXX'R{X) T d. 



7r 



2X 

oo 



I dii?«(i) Hds- 



smUi-'s 



(6.17) 



If desired, the integral over s in (I6.17P can be expressed as a sum of elementary functions 
and the sine integral function. The form in (I6.17P is however more convenient for the 
observations that we shall make below. 

Combining (I6.15P and (16.170 . and using the linearised Einstein equation (16. 6p to 
write = STrGp, where p is the matter density, the gravitational correction to the 
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Im(s) 




Re(s) 



Figure 1: The contour in the complex s plane for evaluating the s-integral in fl6.16p for 
u; < 0. The branch point at s = 2X is indicated by a cross and the cut is at s > 2X. 



Minkowski space transition rate takes the final form 

sin(2cjX) 



dxp(x) 



(2e 



X2 



ds 



smu^s 



2X 



(6.18) 



where we recall that X is defined by X := |x — x|. As a check, we note that the integral 
in fl6.18p converges in absolute value: the quantity in the brackets is is of order 0(X~2) 
as X ^ oo and of order 0(X~^) as X — >• 0, and p(x) = 0(|x|~^) at |x| oo. 

Two observations on the gravitational correction f l6.18p are immediate. First, the 
correction vanishes for u > 0. The excitation rate is thus zero: a static Newtonian 
gravitational field causes no excitations in a static detector. 

Second, the correction depends on the matter distribution within the source even 
when the source is compact and we consider the leading behaviour far from the source, 
|x| — oo. Classical intuition might suggest that in this limit the de-excitation rate 
should only depend on the monopole moment of the mass distribution, the total mass. 
However, as the second term in the brackets in fl6.18p is of order 0(X~'^) at X ^ oo, 
the leading contribution at |x| oo comes entirely from the first term in the brackets 
and can be evaluated, with the result 



G(2e- 1) 
27r|xP 



cos 



sin(2ci;|x|) / dxp(x) cos(2cx;x ■ x) 



+ 0( 



(2c(j|x|) J dxp(x) sin(2c<jx-x) 



') 



(6.19) 



where x := x/|x|. If \ujRo\ <^ 1, where Rq is the characteristic length scale of the 
source, the square brackets in f l6.19p can be approximated by Msin(2u;|x|), where M 
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is the total mass, but outside this hmit the square brackets depend also on the higher 
multipole moments of the source. It is worth mentioning that when Rq is a typical 
stellar or planetary scale and is a typical atomic frequency, we in fact have |a;i?o| ^ I5 
in which limit the integrals in (16.191) can be estimated by WKB techniques. The de- 
excitation rate of an atom far from a star therefore carries an imprint of the internal 
structure of the star. 

6.2 Constant density star 

As an example, we consider a spherical star of constant density po and radius Rq. 
The Ricci scalar has now a discontinuity at the surface of the star, and this spacetime 
therefore falls outside the smooth setting in which we have been working. We suspect 
that the non-smoothness is a technical issue that will not have a significant effect on the 
particle detector and shall proceed, with due caution. 

The integrals in f l6.18p can be evaluated in terms of the sine integral function. When 
X is outside the star, r := |x| > Rq, we find 

AFr{uj) = ^ ^2cu(r + Ro) cos[2tu(r - i?o)] - 2cj(r - i?o) cos[2cu(r + i?o)] 

+ sin[2u;(r - Rq)] - sin[2cL^(r + Rq)] 
+ 4u}\r^ - Rl) (^Si[2cu(r - Rq)] - Si[2cu(r - i?o)]) 

+ ^ ^647ra;Vi?^ - (2u){r + 3i?o) - 4u;^(r^ + r'^R^ - hrRl + 3i?i^)) cos[2cj(r - i?o)] 

+ {2u{r - 3i?o) + 4cu=^(-r3 + r'^R^ + Sri?^ + 3i?3) j cos[2c<;(r + Rq)\ 
+ - 3 + 2cj2(r2 + 2ri?o - 3i?^)) sin[2cj(r - i?o)] 
+ (^3 + 2cj2(-r^ + 2ri?o - 3i?o)) sin[2cj(r + i?o)] 

+ 8cj^(r - Rof{r + 3i?o)Si[2cu(r - i?o)] - 8tu^(r + R^fir - 3i?o)Si[2tu(r + i?o)] 

(6.20) 

Consider in particular the limit in which x approaches the surface of the star, r ^ R^: 
this presumably is the situation with the best experimental prospects of observing the 
gravitational correction to the de-excitation rate. We find 

AFr,{^) = -p^e{-iu) huRo{8uj''Rl - 1) cos(4cui?o) + (3 - 12^ + Slu^RD sin(4cji?o) 
Aoiu Rq L 

+ 8ujRo{87iuj^Rl + 6^ - 1) + 128a;^i?^ Si(2cji?o)l • (6.21) 
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For \ujRo\ ^ 1, the asymptotic behaviour of fl6.2ip is 



AFn,{uj)r.^(^-l]Q{-uj). (6.22) 



\ 6 

If po is the density of the Earth and u is a. typical atomic frequency, the ratio of (16.220 
to the Minkowski vacuum transition rate [to'/(27r)]0(— cu) is of order {Gpo/u'^) ~ 10^^^. 
If the transition rate for the electromagnetic field behaves qualitatively similarly to that 
in our scalar field model, we conclude that the gravitational correction to decay rates in 
atomic physics laboratory experiments is unobservably small. 



7 Conclusions 

In this paper we have discussed the instantaneous transition rate of an Unruh-DeWitt 
detector that is coupled to a scalar field in an arbitrary Hadamard state in curved 
spacetime. We started with a detector that is switched on and off smoothly, in which 
case the detector's response function F^u) is well defined and can be expressed as the 
integral of a Wightman distribution with a standard ie regulator. We showed that the 
limit e ^ can be taken explicitly, with the result 

00 1 f°° ds 

^H = -^ J du[xiu)f + ^ J dux{u)[xiu) ~xiu- s)] 



/OO fOC 
duxiu) / dsxiu — s)Re 
OO Jo 



(7.1) 



where x is the switching function and Wq is the pull-back of the Wightman distribution 
to the detector's world line, with the e-regulator having been taken pointwise to zero. 
We then showed that when the switch-on and switch-off have a fixed shape but take 
each place within the time interval S, the sharp switching limit 6 —>■ results into a 
logarithmic divergence in F{uj), but the derivative of F{uj) with respect to the total 
detection time At remains finite and is given by 



47r 



ds Re 



Wo(r,r-s) + 



+ 



1 



27r2Ar 



+ 0(5). (7.2) 



As a consequence, the difference AFt-{uj) in the transition rates of two detectors in 
different quantum states, on different trajectories and even in different spacetimes, but 
having the same switching function, has a well defined 5 — > limit, given by 



AF,(cu) = 2Re / dse"^-^^ [W^ 



r, r 



0]- 



(7.3) 



The case of a detector switched on in the infinite past can be defined by the Ar — »• cxo 
limit in (17. 2p and (17.31) . subject to suitable asymptotic conditions. We emphasise that 
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all the integrals in the above formulas are integrals of ordinary functions, no longer 
involving ie regulators or other distributional aspects. These results generalise to the 
setting of general Hadamard states in curved spacetime the results that were obtained 
for the massless field in the Minkowski vacuum in f8l. 



We applied the difference formula (17.31) to two situations in which the reference state 
can be conveniently chosen to be an inertial detector in the Minkowski vacuum. First, 
we considered an inertial detector coupled to a massless field in the Rindler vacuum 
in Minkowski space, finding that the transition rate diverges logarithmically as the 
detector approaches the Rindler horizon. Second, we considered a detector at rest in a 
static, Newtonian gravitational field, coupled to a massless field with arbitrary curvature 
coupling, in the Boulware-like vacuum defined with respect to the global timelike Killing 
vector. We found the excitation rate to be zero, but the de-excitation rate acquires a 
gravitational correction that depends on the details of the mass distribution within the 
source, even in the limit in which the source is compact and the detector is far from the 
source. Using a spherical constant density mass distribution as an example, we estimated 
the gravitational corrections to decay rates in atomic physics laboratory experiments on 
the surface of the Earth to be suppressed by 42 orders of magnitude. 

A technical assumption throughout the paper was that both the spacetime and the 
detector trajectory were taken smooth. Given that the final formulas (17.11) and (17.21) 
remain well defined whenever the trajectory is sufficiently differentiable for the term 
to subtract the non-integrable part in Wo, it is tempting to suspect that the smoothness 
assumption on the trajectory could be relaxed. To investigate this question, two steps 
would need to be addressed. First, in section [2] we justified the use of the ie-regulator in 
the pull-back of the Wightman distribution in (12. 6p by the theorems of [T3] , which 
are formulated for smooth submanifolds: how do these theorems generalise to a lower 
degree of differentiabilityllll Second, to obtain in section [3] the estimates in (13.111) . we 
assumed a trajectory that is and has a suitably bounded remainder term in the 
Taylor expansion (say, would suffice): could the techniques of section [H] be improved 
to relax this assumption? 

Our results provide tools for investigating a particle detector's response in time- 
dependent situations in curved spacetime. One set of questions with which these tools 
could prove useful are thermal effects in black hole spacetimes in the time-dependent 
setting. For example, a detector that is falling freely into a static black hole in a 
Boulware-type vacuum [22] would be expected to have a divergent response at the 
horizon, in analogy with the Rindler vacuum analysis in our section [5l but might any 
thermal characteristics survive in the response of a detector falling through the horizon 
in an Unruh-type state [Ij or in a Hartle-Hawking-Israel type state [U |30]? From a 
complementary angle, consider the spacetime of a collapsing star, in a quantum state 
that was Boulware-type in the distant past: how does the response of a detector at a fixed 
position outside the star evolve from that found in section [6] to the thermal response? 
In particular, what is the time scale of this evolution, and how is the detector's response 

^This issue arises already in the Minkowski vacuum analysis in [S]. 
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in this situation related to the the outoing energy flux that develops, or to any notion 
of 'particles' in the associated Bogoliubov transformation? On a more speculative note, 
might there be a relation between the response of a detector and dynamical or evolving 
horizons [31] ? 
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A Appendix: Divergence of the transition rate at 
the Rindler horizon 



In this appendix we verify the logarithmically divergent asymptotic form (15.71) for the 
difference AFt-{uj) of the inertial detector transition rates in the Rindler vacuum and 
the Minkowski vacuum as the detector approaches the Rindler horizon. 

We start from formula (15.61) for AFt-{uj). Writing r = X{1 — e) and scaling the 
integration variable in (15. 6p by s — » Xs, we have 



AFriu) 



1 



27r2X 



"-^ , cos(/?s) 
ds — 



1 + 



1 



+ 



1 



2(l-e)-s lln(l + f) In (1-2^) 



(A.l) 



where f3 := uoX and a := 1 — (tq/X). We shall find the asymptotic form of ( lA.ll) for 
e ^ 0+ with fixed X, a and j3. Note that < a < 2. 

We note first that we may replace the upper limit of the integral in (lA.ip by a 
at the expense of an error of order 0(e). To handle the remaining integral, we split 
the interval (0, a) into the subintervals (0,?;) and (?],«), where rj := ^/e. We denote 
the contributions from the two subintervals by respectively {2'k'^X) Ii and {2'n'^X) I2 
and provide separate estimates for each. A key tool for controlling the logarithms will 
be the Laurent expansion 



ln(l +x) 



(A.2) 



Consider Ji. We introduce the new integration variable r := s/i], with the range 
< r < 1. We replace the second logarithm in the integrand by the first two terms 
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m f[0) . at the expense of an error of order 0{rj) in /i, and find 



Jo r [2{1 - T]^) - r]r\ \^inh + i j r 

As the quantity in the large parentheses in flA.Sp is bounded, we may make the replace- 
ments cos{f3T]r) — i> 1 and [2(1 — rj'^) — ?7r] ^ — > [2(1 — 77^)] ^ at the expense of respective 
errors of order 0{rf) and 0{ri). Changing the integration variable to y := r/rj, we then 
obtain 

We now concentrate on the divergent part. With errors of order 0(1), we first replace 
the lower limit of integration in flA.4p by 1, then drop the |/~^ term and make the re- 
placement [yln (1 + y)]^^ [(1 + y) In (1 + y)]~'^- The remaining integral is elementary, 
with the result 

/i = |lne + iln(-lne) +0(1). (A.5) 

Consider then I2. In the integrand shown in (lA.ip . we replace the argument of the 
second logarithm by (l — f ) , the fraction s/ [2(1 — e) — s] by s/ (2 — s) and the argument 
of the first logarithm by s/e, with elementary estimates showing that each step produces 
in I2 an error of order 0{rj), and obtain 



° , cos(/?s) 
as 



n 



S2 



+ 



2-s \\n{s/ri^) In (l - f ) 



+ 0(r/). (A.6) 



The terms in (I A. 61) that do not involve In^s/r^^) can be handled by a straightforward 
small s Laurent expansion, and their contribution to I2 is jlnr] + 0(1). In the term 
involving \ia(^s / 1]'^^ , the replacements cos(/5s) — > 1 and s/(2 — s) ^ s/2 can be verified 
to produce in I2 errors of order 0(l/(ln?7)), and the remaining integral is elementary 
and evaluates to 0(1). Hence 

J2 = |lne + 0(l). (A.7) 
Combining these results, we have 

AF.(o.) = ^[ilne + iln(-lne) + 0(l)] . (A.8) 
Formula fl5.7p follows by substituting e = 1 — (r/X). 
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